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J> ' Abstract 

in 

This paper discusses a new model for galactic dark matter by combining an anisotropic 
pressure field corresponding to normal matter and a quintessence dark energy field 
having a characteristic parameter Ug such that — 1 < Ug < — ^ . Stable stellar orbits 
together with an attractive gravity exist only if ojq is extremely close to — ^ , a result 
consistent with the special case studied by Guzman et al. (2003). Less exceptional 
forms of quintessence dark energy do not yield the desired stable orbits and are 
therefore unsuitable for modeling dark matter. 

X' 

- ■ 1 Intoduction 

The existence of dark matter has been suspected since the 1930's due to stellar motions in 
the outer regions of galaxies [H El E] . Zwicky [U [5] discovered the presence of dark matter 
on a much larger scale by studying galactic clusters. The presence of such nonluminous 
matter was subsequently confirmed by observing the flatness of galactic rotation curves 
[TJ m [9l [TOl [11]. While the dark matter provides the necessary gravitational field to 
account for the observation, its composition remains unknown, probably consisting of 
particles not included in the standard model. 
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The origin of the dark-matter problem is the measurement of the tangential velocity v''' 
of stable circular orbits of hydrogen clouds in the outer regions of the halo. The observed 
constant velocity can only be explained by assuming that the energy density decreases 
with the distance as so that the increase in the mass of galaxies is proportional to the 
distance r from the center. 

Of the many proposed candidates for dark matter, the most favored is the standard 
cold dark matter (SCDM) paradigm [121 [13]. Another favorite is the A CDM model that 
is related to the accelerated expansion of the Universe [Til [15] . For a summary of some 
alternative theories such as scalar-tensor and brane-world models, see Rahaman et al. 
[mini [18]. 

It is now recognized that on a cosmological scale, our Universe is pervaded by a 
dynamic dark energy that causes the expansion of the Universe to accelerate: a/ a = 
— ^y^(p -|- 3p). In the equation of state p = wp, the range of values —l<w< —1/3, 
resulting in a > 0, is usually referred to as quintessence dark energy. Moving to the 
galactic level, however, there is no model consistent with the cosmological one to help 
explain the nature of dark matter. 

In this paper we propose a new model based on a combination of a quintessence-like 
field and a field with (possibly) anisotropic pressure, representing normal matter. We 
study conditions leading to stable circular orbits and attractive gravity, as well as the 
observed equation of state of the galactic halo, based on a combination of rotation curves 
and lensing measurements [19]. The conditions required for stable orbits and attractive 
gravity are rather stringent, but by comparing global and local effects, they proved to be 
consistent with earlier results obtained by Guzman et al. [20] . 



2 The model 

In this paper the metric for a static spherically symmetric spacetime is taken as 

ds^ = -e^'^'^dt^ + e^^'Ur^ + r^dO^ + sin^O dcf)''), (1) 

where the functions of the radial coordinate r, i/(r) and A(r), are the metric potentials. 

Now we consider a model which contains a quintessence field and a second field with 
an anisotropic pressure representing normal matter. Here the Einstein equations can be 
written as 

G^, = 87rG(T^, + r^,), (2) 

where r^,^ is the energy momentum tensor of the quintessence-like field, which is character- 
ized by a free parameter Ug with the restriction — 1 < < — |. According to Kiselev [21] . 
the components of this tensor need to satisfy the conditions of additivity and linearity. 
Taking into account the different signatures used in the line elements, the components 
can be stated as follows: 

r't = <■ = -P„ (3) 
r'e=4 = l{3u;, + l)pg. (4) 
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Moreover, the most general energy momentum tensor compatible with spherically sym- 
metry is 

Ti: = (p + Pt)u''u, - ptgii + {pr - ptKVu (5) 
with u^^u^ = — 1. For our metric the Einstein field equations are stated next. 
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87rG'(p + p,), 



— = S-nGiyPr - pq) 



"1 z/'" 

- + - 

l(z.r + ^"-^AV + l(z.'-A') 



(6) 
(7) 



(8) 



3 The solutions 

It is well known [TTj |22] that the fiat rotation condition gives the solution 

= Bor\ 



(9) 



where / = 2v'^'^, v"^ is the rotational velocity, and Bq is an integration constant. According 
to Ref. |23], the observed rotational curve profile in the dark matter dominated region is 
such that the rotational velocity v'^ becomes more or less constant with f*^ ~ 300 km/s 
(~ 10~^) for a typical galaxy. So I = 0.000001, as in Nandi et al. ^1]. We also assume 
fairly large distances in Kpc. (Of course, certain galaxies, particularly disk and dwarf 
galaxies, are not typical in the above sense. For further discussion, see Refs. [25| 126].) 
To continue the analysis, we use the following equation of state (EoS): 



Pr 



mp, m > 0, 



(10) 



where m is a parameter corresponding to normal matter. In the case of an anisotropic 
pressure, the radial and lateral pressures are no longer equal. So let us assume the form 



Pt = ap, 



a > 0, 



:iii 



for the lateral pressure. In this manner we obtain simple linear relationships between 
pressure and energy- density, but with pr not equal to pt, unless, of course, m = a. 

Now, from equations (l6l) to (IE]) and then using equations (l9l), (fTOj) . and (fTTj) . one 
obtains the following simplified form (after some manipulation): 



+ 



ae 



(12) 
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m 


a 


a 


c 


-0.33333378 


0.1 


0.1 


-1.287692503 x lO'^ 


-1.133845791 x 10" 


-6 




0.1 


0.3 


-1.220000166 X 10-6 


-8.670586114 x 10" 


-7 




0.1 


0.5 


-1.178095381 X 10"^ 


-7.019046229 x 10" 


-7 




0.3 


0.5 


-1.192174084 X 10"^ 


-7.573912227 x 10" 


-7 




0.5 


0.5 


-1.204000197 X 10"6 


-8.039999743 x 10" 


-7 




0.5 


0.8 


-1.161613061 X 10-*^ 


-6.454838537 x 10" 


~7 




0.7 


0.7 


-1.183161478 X 10-*^ 


-7.315484069 x 10^ 


-7 




0.8 


0.8 


-1.176823715 x 10"^ 


-7.062353288 x 10" 


-7 


-0.4 


0.1 


0.1 


-0.1718752386 


-0.1718749261 






0.1 


0.3 


-0.1309528143 


-0.1309523381 






0.5 


0.5 


-0.1250004609 


-0.1249999609 




-0.5 


0.1 


0.1 


-0.4400003664 


-0.4399998064 






0.1 


0.3 


-0.3333338889 


-0.3333332222 






0.5 


0.5 


-0.3333338889 


-0.3333332222 




-0.8 


0.1 


0.1 


-1.327586957 


-1.327585577 






0.1 


0.3 


-0.9871804212 


-0.9871791391 






0.5 


0.5 


-1.166667514 


-1.166666181 





Table 1: Table showing typical values of a and c corresponding to various values of oOg, 
m, and a, with special emphasis on Ug = —0.33333378. 



where 



and 



Ul + l){3ujg + 1) + f - l[{3ujg + l)m + 2a]/[2(l + m)] 



Eq. (12) yields 



I + I + [i3cOg + l)m + 2a]/[2(l + m)] 

|(3a;, + i) 

I + i + [{3ujg + l)m + 2a]/[2(l + m)] ' 



c D 
- + — . 

a 



(13) 



(14) 



(15) 



where D is an integration constant. This constant would normally be obtained from the 
junction conditions, but the exact boundary of the galactic halo is unknown. In some 
situations D may have to be restricted for physical reasons. 

Both a and c depend on Uq, as well as on m and a. To get an overview, some typical 
values are listed in Table 1, using / = 0.000001. (The special case u = —0.33333378 will 
be discussed in the next section.) 

From Eqs. ([7]) and ([8]), one can readily obtain the following solutions for p, pg-. 



P 



1 

87rG(l + m) 



'm 



r.a+2 



+ 



Ic 



(16) 



4 



and 



8nG 



[aD ~ D - {ID + Da)/{l + m)] 1/c 



Ic 



SttG 



D(a-l) ^ 1 



a(l + m) 
D(a + /) 



ra+2 



To complete the analysis below, we also need the following additional results: 



p{effective) = p + pq 



8nG 



[Gt 



SttG 



D(a 



r 



a+2 



(17) 



:i8) 



Pr{ef fective) = Pr - Pq 



SnG 



[Gr 



8nG 



+ fc D\ 1 



(19) 



Pt{ef fective) 



(3avfl) 
+ ^ Pq 



8nG 



8nG 
c 



[Ge 



+ 



D 



2ar2 2r'^+2 



Da'' f l_ 
~2 ~ or'' + Da \ ^2 



(20) 



and 



SttG [p{ef fective) + pr{ef fective) + 2pt{ef fective)] 
^Dia-l ' ' 



r.a+2 



1 

+ ^ 



D fc D\ /, P .1 , 



(21) 



4 Stability of circular orbits 



Let t/" = ^ be the four-velocity of a test particle moving solely in the space of the 
halo and suppose we restrict ourselves to 9 = 7r/2. Then following Nandi et al. [21], the 
equation g^^U'^U'^ = —rn^ can be cast in a Newtonian form 



which gives 

V{r) = 
Here the constants 



{c/a + D/r") 



(22) 





"c D' 




+ 


- + — 


(-^)} 




a r" 





E = and L = 

mo mo 



(23) 



(24) 
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are, respectively, the conserved relativistic energy and angular momentum per unit rest 
mass of the test particle. For circular orbits we have r = R, a constant, so that ^ = 0. 

=ij= 0. From these two conditions follow the conserved parameters L and E: 



Also, f 



L = ± 



2-1 



R- 



(25) 



using L in V{R) = —E"^, we get 



E = ± 



(26) 



The orbits will be stable if \r=R< and unstable if \r=R> 0. Putting the ex- 
pressions for L and E in \r=R, we obtain, after straightforward calculations, the final 
result: 



c 21 



r=R 



a2- I 



n 1^ 1 3 



+ 



2 

2^1 
21 



R\ 
a(a + 1)D 



1 



1 



:D 



I 



il + l] 



1 3 



21 



+ 



2-1 



{2a)D 



^a+l+l Ra+2 



(27) 



To allow a comparison to p [Eq. f|T6|) ]. V"{R) has to be rewritten. Before doing so, we 
recall that a and c are both negative, while c/a is less than 1. Furthermore, a and c 
depend on Ug and can be made small by choosing Ug close to —1/3. But \a\ must be 
greater than Z; otherwise the negative second row in Eq. (27) cannot catch up with the 
positive first row. Of course, \a\ must be less than 1, so that a{a + 1) is negative. 



V"{R) 



2 c 



I R a 



il + l] 



Ri 



+ 



+ 



3 
R 
2 



2-lR 
1 2 
2-lR^ 



Daia + 11 



-Dl 



{i + r 



' 1 

1 



+ 



1 

'r 

3" 

'r 

1 



1 
1 
2 



2-lR 



Di2a)l 



1 

'R} 



1 

'r 



1 



(28) 



For quintessence dark energy, -I <u)q< -\. The case where ujg is extremely close 
to — I will be referred to as marginal quintessence dark energy, otherwise as ordinary 
quintessence dark energy. 
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4.1 Ordinary quintessence dark energy 

Suppose uj is not real close to — |, so that \a\ is not real close to / and a + / is not real 
close to zero (Table 1). Then in Eq. ( 128|) . the second term in each bracket is relatively 
small and we can concentrate on the first term in each case: 



V"{R) 



1 2 1 
2^1 RE} 



c 
a 

1 2 1 



l{l + 1) + [Da{a + 1) + Dl{l + 1] 

' D 



2alD] 



1 



2 - / i? i?' 



y + a + r + l + 2al) + + 1) 



1 2 1 



2-lRR^ 

Now compare this equation to Eq. (16) for r = R: 

- 1 1 \D_ c 



^Ja + lf + {a + l)] + -J{l + l) 



(29) 



(30) 



Since / = 0.000001, Z(Z + 1) ~ /. To make physical sense, D has to be small enough to 
make p positive. But then V"{R) is a fortiori positive because of the extra positive term 
(a + ly. In other words, we do not have a stable orbit. Conversely, if the orbit is stable, 
then p has to be negative. 



4.2 Marginal quintessence dark energy - global and local effects 

As already noted, referring to Table 1 and Eq. (fT3l) . as Uq gets close to — |, \a\ gets close 
to / and a + / close to zero. Then the second term inside each bracket in Eq. fl28p is no 
longer negligible. Collecting these terms, we get 



1 2 
2-lR 





±1 







(31) 



Observe next that both and al are negligible. Supppose that in Eq. ( 129|) . D\a+l\ < (^)Z; 
(a + 1)^ is negligible and is close to unity for R between 100 and 500 Kpc. So in Eq. 
(l3T|l . since D{—a — I) < the quantity 



-/(-3) + L>(-a - / - 2/) 

a 

is less than — 2(-)/ — 2DI and hence negative. In other words, we now have p > and 
V"{R) < 0, finally yielding a stable orbit. 

Remark: A simple alternative argument is the following: consider the limiting case 
Uq — — |. In that case, c — )■ and \a\ < I [from Eqs. ( [T3|) and (fT4l) ] and the conclusion 
follows from Eqs. (150]) and (ET]). 

An example of the marginal quintessence-energy case is discussed in Guzman et al. 
[20] . It is assumed that the luminous matter does not contribute significantly to the 
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energy density of the halo and may be disregarded, leaving an isotropic perfect fluid. It 
is shown that the equation of state may be taken to be 

1+^20 

or, equivalently, ujq ~ —0.33333378, which is indeed close to — |. Table 1 shows that a is 
small enough to meet the requirements in Eqs. fl30l) and fl3Tl) for both the isotropic and 
anisotropic cases for R between 100 and 500 Kpc. 

While one usually thinks of quintessence dark energy as a cosmological phenomenon, 
in particular, as the cause of the accelerated expansion, the above results, taken together, 
point to a rather different effect on the galactic level. The marginal quintessence case, 
being just over the line into quintessence territory, yields stable orbits in our combined 
model. Ordinarily quintessence dark energy is unable to do so: evidently, its repulsive 
action is too strong. The difference in behavior is reminiscent of the expansion itself: 
while the Universe expands on a cosmological scale, bound systems, such as galaxies, are 
completely unaffected. In other words, the same phenomenon has different effects on the 
cosmological and local levels. 

The idea that dark matter and dark energy are connected is not new. In discussing 
Chaplygin gas, it is argued in Ref. that dark matter and dark energy ought to be 
different manifestations of the same entity. 

Returning to some of the other conditions in Sec. [3l in the marginal quintessence 
case, both pq and p{ef fective) are positive as long as D is small enough, say 0.1, while 

= —pq < 0, as required. For the same D, pr{ef fective) is likely to be negative 
(depending on the precise values of D, a, and c), while pt{ef fective) is positive for any 
D. 



5 Attraction and total gravitational energy 

Suppose we now consider the question of attractive gravity by examining the geodesic 
equation for a test particle that is moving along a circular path of radius r = R: 



dr^ 



M7 



dr dr 



0. 



(32) 



This equation implies that 



d'^r 
d^ 



D 



Da 



D 



f dr 



\dT 



(33) 



As before, as long as ^ = 0, the expression is negative and we conclude that objects are 
attracted toward the center. But this would only be relevant in the marginal quintessence 
case. For this case we can also determine the total gravitational energy Eg between two 
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fixed radii, ri and r2 [28]: 



Eg = M- Em 



An 





rr2 












ri 












1 




I T\ 







[1 - Ve^]pr'^dr 



Here 



c/a + D/r^ 



M = A-K 



8nG 



D{a-1) (1-c/a) 



r.a+2 



r2 



pr'^dr 



ri 



is the Newtonian mass given by 



M = Att 



r2 



pr'^dr 



G 



c/a)r 



D 



2r 



a-l 



r2 



J ri 



r^dr. (34) 



(35) 



(36) 



Em is the sum of the other forms of energy such as the rest energy, kinetic energy, and 
internal energy. So the total gravitational energy is given by 



E. 



1 

G 



c/a)r 



D 



2r 



a-l 



c^rF[(0.5, 1/a); (1 + 1/a); -ar^D/c] 



r2 



ri 



1 



Dd - ^y(-o.5a+i) + l/g A5); (0.5 + 1/a); -crV(aD)] 

Da(-0.5 + 1/a) 



^■2 



• (37) 



ri 



Letting D = 0.1 again, let us plot the total gravitational energy, assuming a distant 
halo region, typically r ~ 200 Kpc, shown in Fig. 1. The value of Eg is small but negative, 
showing that gravity in the halo is indeed attractive for the marginal quintessence case. 
(As before, this is not the case for ordinary quintessence dark energy.) 



6 The observed equation of state 

For the marginal quintessence case, the equation of state of the halo fluid can be obtained 
from a combination of rotation curves and lensing measurements. But first we need to 
rewrite the metric, Eq. ([T]), in the following form: 



ds^ = -e^^i^^dt' + 



1 — 2m{r)/r 



dr^ + r'^dnl, 



where 



and 



- [ln5o + /lnr] 



(38) 



(39) 



m{r) 



c 
a 



D 

„a—l 



(40) 
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Figure 1: The variation of Eg with r in Kpc. The lower hmit of integration in Eq. ( !24l) is 
fixed at ri = 100 Kpc, while r2 varies from 100 to 500 Kpc. We choose G = 1, D = 0.1, 
and ~ lO^'^ (300km/s) for a typical galaxy. 



As discussed in Ref. [21], the functions are determined indirectly from certain lensing 
measurements defined by 



lens 



$fr) 1 f mir) , lni?o / , 1, c , D 

-^;—dr = — |--lnr + -lnr — — lnr + 



and 



miens = ^r2$'(r) + ^m(r) = j + ^ 



4a 



cr D 



Aar°- 



a r 



Of particular interest to us is the dimensionless quantity 

due to Faber and Visser [19]. The subscript RC refers to the rotation curve; thus 



(41) 
(42) 

(43) 



W = 0(r) = -[ln5o + /lnr] 



and 



niRc = r $'(r) 



Ir 



(44) 
(45) 



The primes denote the derivatives with respect to r. The resulting expression is 
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m 



lens 



1 I {al — a + c)r"- — Da{a — 1) 
3 '^mlns - ^'rc 3 [ Da{a - 1) + (a - c)r'^ 

1 Da{a 

3 Da{a — 1 



1) + (a — c)r'^ — air"- 



(46) 



We conclude that for our marginal quintessence case, uj{r) is only minutely less than — |, 
which is consistent with our earlier result. 



7 Conclusion 

This paper discusses a new model for galactic dark matter by combining a quintessence 
field with a normal matter field that may have an anisotropic pressure. The quintessence 
field is characterized by a free parameter Uq with the restriction — 1 < < — |. For con- 
venience, the case where Uq is extremely close to — | is referred to as marginal quintessence 
dark energy, otherwise ordinary quintessence dark energy. It was found that for ordinary 
quintessence, stable orbits exists only if p < 0. For the marginal quintessence case, how- 
ever, stable orbits do exist for p > 0, together with an attractive gravity. This result 
is consistent with that of Guzman et al. [20], who employed a perfect-fluid model for 
dark matter, while disregarding ordinary matter. The equation of state of the halo field, 
obtained from a combination of rotation curves and lensing measurements, is also con- 
sistent with the marginal case. In summary, quintessence dark energy appears to have 
different effects on the cosmological and local levels: while quintessence dark energy is 
responsible for the accelerated expansion, marginal quintessence dark energy is suitable 
for modeling dark matter, although ordinary quintessence dark energy is not. Finally, it 
has been argued that dark matter and dark energy ought to be different manifestations 
of the same entity. 
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